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STM microscopy of the CDW in lT-TiSe2 in the presence of single atom defects 
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We present a detailed low temperature scanning tunneling microscopy study of the commensurate 
charge density wave (CDW) in ir-TiSe 2 in the presence of single atom defects. We find no significant 
modification of the CDW lattice in single crystals with native defects concentrations where some bulk 
probes already measure substantial reductions in the CDW phase transition signature. Systematic 
analysis of STM micrographs combined with density functional theory modelling of atomic defect 
patterns indicate that the observed CDW modulation lies in the Se surface layer. The defect patterns 
clearly show there are no 2JT-polytype inclusions in the CDW phase, as previously found at room 
temperature [Titov A.N. et al, Phys. Sol. State 53, 1073 (2011). They further provide an alternative 
explanation for the chiral Friedel oscillations recently reported in this compound [J. Ishioka et al., 

Phys. Rev. B 84, 245125, (2011)]. 
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The transition metal dichalcogenide (TMD) lT-TiSe 2 
has kept the scientific community wondering about a 
number of its striking physical properties for more than 
four decades [ll-0|- ir-TiSe 2 is a layered compound con¬ 
sisting of a hexagonal layer of Ti sandwiched between 
two hexagonal layers of Se to form Se-Ti-Se sandwiches 
that stack via weak Van-der-Waals (VdW) forces to form 
a single crystal. The band structure of lT-TiSe 2 , as de¬ 
termined by angle-resolved photoemission spectroscopy, 
consists primarily of a Se 4p-valence band at the T point 
and a Ti 3d-conduction band at the L point of the Bril- 
louin zone. But it is still debated whether it is a semi¬ 
conductor or a semimetal with evidences claimed for both 
alternatives 

Below Tqdw ~ 202Ar, ir-TiSe 2 undergoes a second 
order phase transition into a commensurate charge den¬ 
sity wave (CDW). A comprehensive theory of this CDW 
formation is yet to be developed. Two main mecha¬ 
nisms are currently considered, driven either by a Jahn- 
Teller distortion E El or an excitonic ground state 
Eilillll. The CDW phase has been found to melt 
mon copper intercalation Q or when applying pressure 
[7|. In both instances, superconductivity develops in a 
dome shaped region around some optimal doping or op¬ 
timal pressure, with a maximum critical temperature of 
4.IK and 1.8K, respectively. More recently, chiral prop¬ 
erties have been reported for the CDW in pristine and 
copper intercalated lT-TiSe 2 based on polarized optical 
reflectometry and scanning tunneling microscopy (STM) 

[IMl. 

Here we focus on the CDW instability in ir-TiSe 2 in 
the presence of native atomic scale defects. Past studies 


performed using macroscopic probes including resistiv¬ 
ity, magnetic susceptibility and optical reflectivity have 
found atomic intercalation and substitution to be detri¬ 
mental to the CDW E1I3- This compound is usually 
non-stoichiometric with a strong correlation between in¬ 
creasing crystal growth temperature and Ti self-doping 
leading to the collapse of the CDW phase transition 
signature in temperature dependent resistivity measure¬ 
ments [H . STM offers new opportunities in allowing the 
simultaneous mapping of individual single atom defects 
and the CDW in real space, as well as measuring the lo¬ 
cal density of states (LDOS) around the Fermi level by 
tunneling spectroscopy. This technique has revealed a 
distorted CDW superlattice in doped ir-TaS 2 [sm, 
a one-dimensional CDW in calcium intercalated graphite 
[ 2 ^ and a finite CDW amplitude in the vicinity of in¬ 
trinsic defects in 2i/-NbSe2 well above the bulk Tcdw 
These examples highlight the possibility to gain in¬ 
sight into the CDW phase and its formation mechanism 
by means of STM in the presence of atomic defects and 
impurities. 

ir-TiSe 2 single crystals were grown by iodine vapour 
transport and cleaved in-situ below 10~^mbar at room 
temperature. All measurements were performed on crys¬ 
tals grown at 650°C, except for the micrograph in fig. 3j 
acquired on a crystal grown at 575°C to better observe 
atomic features unrelated to intercalated Ti. Constant 
current STM micrographs were recorded at 4.7K using 
an Omicron LT-STM and a SPECS JT-STM, with bias 
voltage Vbias applied to the sample. Base pressure was 
in both cases better than 5xl0“^^mbar. DFT model cal¬ 
culations were performed using the plane wave pseudo- 
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FIG. 2. (color online), (a) Ball-and-stick model of the IT- 
TiSe 2 lattice showing the positions of the native defects A-D. 
(b) Model representation of the inequivalent 1/4 (filled sym¬ 
bols) and 3/4 (empty symbols) lattice sites for defects A-D in 
the commensurate CDW phase. 


FIG. 1. (color online). Simultaneously measured empty- 
states (a: Vbias =0.15V) and filled-states (b: Vbias = —0.15V) 
STM micrographs of a lT-TiSe 2 single crystal grown at 650°G. 
Image size: 22.2 nm xll.4 nm, It = 0.2nA, T = 4.7K. Native 
defects are labelled A, B, G, D. 


potential code VASP [22|, l23|, version 5.3.3. Projector- 
augmented waves 241 were used with the PBE [25| ex¬ 
change correlation functional and plane wave cut-offs of 
211eV (ir-TiSe2, I substitutional) and 400eV (O). The 
cell size of our model was 28.035Ax28.035A. The IT- 
TiSe 2 surface was modeled with two layers and the bot¬ 
tom Se layer fixed. A Monkhorst-Pack mesh with 2x2x1 
/c-points was used to sample the Brillouin zone of the 
cell. The parameters gave energy difference convergence 
better than 0.01 eV. During structural relaxations, a tol¬ 
erance of 0.03eV/A was applied. STM ima ges were gen¬ 
erated using the Tersoff-Hamann approach |2fil | in which 
the current I{V) measured in STM is proportional to the 
integrated LDOS of the surface using the bSKAN code 


27|. 


Fig. [T] shows two high resolution STM micrographs of 
lT-TiSe 2 obtained at T = 4.7K with exactly the same 
tip at positive and negative sample bias [ 2 ^. The bias 
voltages of ±150mV have been chosen to enable the si¬ 
multaneous resolution of the 2ao x 2bo CDW reconstruc¬ 
tion on the selenium layer and atomic lattice features 
at opposite polarities. Defects (A-D) correspond to the 
dominant native atomic defects in ir-TiSe 2 identified in 
a recent STM/DFT study based on images recorded at a 
larger bias voltage where the CDW is not resolved (29| . 
These defects are Se surface vacancies (A), iodine (B) 
and oxygen (C) substitution for bulk Se and titanium in¬ 
tercalated into the VdW gap (D). Their positions in the 


lattice unit-cell are shown in Fig. Ha). 

lT-TiSe 2 cleaves between the weakly VdW bonded Se- 
Ti-Se sandwiches, thus exposing a hexagonal Se layer to 
the surface. DFT modeling enables us to identify the 
atomic lattice seen in STM maps with the Se surface layer 
by assigning the observed vacancies (defect A) to missing 
Se surface atoms [ 2 ^. Thus, the commensurate in-plane 
2aox2bo modulation [Fig. [1] is in perfect registry with 
the Se atomic lattice, indicating that the CDW charge 
modulation detected by STM resides in the Se layer. 

As a consequence of the CDW modulation, there are 
two inequivalent sites in the unit-cell for each defect, with 
three times more 2/4 than 1/4 sites [Fig. Hb)]- A survey 
of Se, O and Ti defects in a large area map (50x50nm^, 
23’000 unit cells, 177 defects in total) yields approxi¬ 
mately three times more 3/4 than 1/4 configurations for 
each of them. This uniform statistical distribution of all 
native defects among 3/4 and 1/4 sites implies they do 
not interact strongly with the CDW in this crystal, even 
though its resistive CDW transition is reduced by over 
30% compared to a sample with optimal stoichiometry. 
If they were interacting, we would expect dislocations to 
enable the CDW lattice to accommodate the random de¬ 
fect landscape. Indeed, we find no systematic domain 
formation, dislocations or weakening of the CDW lattice 
due to native defects. This rigidity of the CDW can be 
directly linked to its commensurate nature in ir-TiSe 2 
[^ . In the same survey, we count about 80 intercalated 
Ti (defect D), corresponding to 0.35% self doping, in ex¬ 
cellent agreement with literature for samples grown at 
650°C 0. 

Se vacancies appear as well resolved dark sites inde¬ 
pendent on bias voltage and position at the surface [Figs. 
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FIG. 3. High resolution STM micrographs centered on Se vacancies (a,b), iodine (d,e) and oxygen (g,h) substitutions at 3/4 
sites (a,d,g) and 1/4 sites (b,e,h) with corresponding DFT simulations calculated without including the CDW (c,f,i). (j) Linear 
features around defects A and B observed by STM on lT-TiSe 2 grown (ll.Sx 11.5nm^, Vbias = +150mV, It = 0.2nA). All 
images are taken on crystals grown at 650°C except for image (j) recorded on a sample grown at 575°C to have a reduced 
number of intercalated Ti. 


|3Ua)-(c)]. In contrast, defects B-D are mostly bright and 
best resolved and differentiated at positive Vbias [Fig . [T]. 
Their characteristic patterns revealed by STM [29j are 
slightly modified in the presence of the CDW and de¬ 
pend on their 1/4 or 3/4 configuration (Figs. [3]and|4l). 
Of all defects, iodine substitution for Sedown (defect B) 
is the most difficult to identify. On the 1/4 site, it ap¬ 
pears as a faint enhancement of the three nearest CDW 
maxima [Fig. EJe)]. The 3/4 configuration appears as 
a few atoms long brighter chain extending along one of 
the crystallographic direction [Fig. EKd)]. Similar bright 
atomic chain features are found around defect A [Figs. 
Ha) and (j)]. DFT simulations of defects A and B are 
in good agreement with these experimental observations 
[Figs. He) and (f)]. The linear features are reproduced in 
the model without including the CDW instability. This 
shows they are not a different CDW ground state (e.g. 
ID-CDW), but reflect local strain due to the Se vacancy 
and the larger atomic radius of iodine compared to Se. 
In regions with higher defect densities, these chains co¬ 
operate to form stripy patches in the STM topography, 
but without disrupting the long range coherent CDW 
[Fig. Hi)]- The defects in these regions will produce an 
anisotropic deformation landscape explaining why these 
stripes do not always appear in all three high symmetry 
directions with the same intensity. When the defect den¬ 
sity is low and in the vicinity of intercalated Ti (defect 
D), these stripes are usually less or not visible [Fig. [TJS')]- 
Defects C and D show more complex triangular pat¬ 
terns without the linear atomic features found around 
defects A and B. Oxygen substitution for Sedown (de¬ 
fect C) is characterized by three bright central atoms 
centered on a larger, 60° rotated triangle of three dark 
atoms [Figs. Hs) (h)], in perfect agreement with 
DFT modeling [Fig. Hi)]- Titanium interstitials (defect 


D) appear as two concentric bright triangles centered on 
the defect [Fig. 0] (2^. The central triangles point in 
opposite directions in defect C compared to defect D. 
The triangular outline of defects C and D always point 
in the same direction in a given experiment [Fig. [T], at¬ 
testing the perfect crystalline structure of our ir-TiSe 2 
specimen. The unique triangle orientation and the per¬ 
fect match between the data and the DFT models, which 
were all calculated in the IT-polytype structure, imply 
there are no 2iL-polytype inclusions where the coordi¬ 
nation of the Ti atom changes from octahedral (IT) to 
trigonal-prismatic (2iJ) as found by Titov et al. at room 
temperature 31|. Although this finding cannot exclude 
a Jahn-Teller mechanism for the CDW origin aim , the 
capability to locally identify the IT-phase may become 
instrumental in clarifying the role of local lattice modifi¬ 
cations in the CDW formation. 

The appearance of all native defects, except surface Se 
vacancies (defect A), change slightly depending on their 
1/4 or 3/4 configuration. At positive sample bias, the oxy¬ 
gen substitution (defect C) in the 1/4 configuration to¬ 
tally obscures the amplitude of the three nearest CDW 
maxima [Fig. H^^)] whereas iodine on the same location 
(defect B) enhances them slightly [Fig. H®)]- Interca¬ 
lated titanium (defect D) has an unmistakable triangular 
signature, with very sharp vertices in the 1/4 configura¬ 
tion that become nearly extinct in the 3/4 configuration. 
These different appearances of native atomic defects de¬ 
pending on their configuration (1/4 or 3/4) suggest another 
explanation for the recently reported chiral Friedel oscil¬ 
lations. Our data and DFT modelling show the distinct 
left and right handed patterns discussed by Ishioka et al. 
[l5| correspond in fact to different native defects (O and 
I-substitutions) in the two distinct 1/4 and 3/4 configura¬ 
tions, unrelated to chirality. 
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FIG. 4. High resolution STM micrographs centered on inter¬ 
calated Ti at 3/4 sites (a,e) and 1/4 sites (c,g) with correspond¬ 
ing DFT simulations in the presence of the CDW (b,d,f,h). 
Vbias = -150mV (a,b,c,d) and Vbias = -|-150mV (e,f,g,hi 
It = 0.2nA. 

The native defects are poorly resolved in the negative 
low bias STM micrographs discussed here, except for in¬ 
tercalated Ti (defect D) in the 1/4 configuration and Se 
vacancies (defect A). The dark sites associated with Se 
vacancies (defect A) correspond to holes in the topogra¬ 
phy and are seen as such at both polarities. The other 
defect patterns are primarily electronic and their bias 
polarity dependent visibility observed here is consistent 
with a CDW gap that is biased towards occupied states at 
the Fermi level [16[ . A striking exception to this behavior 
is defect D which is nicely resolved in the 1/4 configura¬ 
tion at Vbias < 0 [Figs. [T][b) and|4l[c)], closely matching 
the DFT modeling. The donor nature of intercalated Ti 
contributing electron states just above the occupied edge 
of the CDW gap [ 2 ^ can explain the finite contrast of 
defect D at negative bias inside the CDW gap. However, 
it is presently not clear why only the 1/4 configuration 
is resolved at Vbias = —150mV [Fig. SKc)] while the ^/i 
configuration remains invisible [Fig. |H[a)]. This ques¬ 
tion demands further investigations, in particular in the 
context of a proposed excitonic ground state Hi- 

In summary, the careful analysis and comparison with 
DFT modeling allows us to assign the surface patterns 
observed in STM micrographs of the CDW phase in IT- 
TiSe 2 exclusively to Se vacancies, O or I-substitutions 
and Ti intercalation [ 2 ^. We have shown the great po¬ 
tential of high resolution STM imaging of the CDW in 
the presence of atomic defects to gain insight into this or¬ 
dered phase. We find that native defects have essentially 
no impact on the CDW lattice, at least up to the level of 
Ti self-doping considered here, where the corresponding 
phase transition is significantly reduced in transport mea¬ 
surements [ij. The only change, besides the characteris¬ 
tic signatures of each defect, is a locally enhanced bright¬ 
ness linked to their doping nature [29j. The theoreti¬ 
cal proposal by McMillan [30[ that atomic defects might 
trigger an incommensurate CDW is clearly not observed 


here. For the first time, comparison with DFT model¬ 
ing further allows us to unambiguously demonstrate that 
the observed defect patterns in the CDW phase are all 
consistent with the IT-polytype, excluding 2iJ-polytype 
inclusions Finally, our study sheds a different light 
on recently published work on the chiral nature of Friedel 
oscillations in the vicinity of defects [1^. We find com¬ 
pelling evidence that the left and right-hand patterns 
identified in that work are the signatures of different IT- 
TiSe 2 native defects located on inequivalent lattice sites 
with respect to the CDW modulation. 
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